This figure "fig 1-1.png" is available in "png" format from: 


http://arXiv.org/ps/astro-ph/9502035vl 



This figure "fig2-l.png" is available in "png" format from: 


http://arXiv.org/ps/astro-ph/9502035vl 





This figure "fig 1-2.png" is available in "png" format from: 


http://arXiv.org/ps/astro-ph/9502035vl 

















This figure "fig 1-3.png" is available in "png" format from: 


http://arXiv.org/ps/astro-ph/9502035vl 



op p p 

p o bi O CJi O QI 



CJi O CJI o 

^ o o o o 

P o o o o o 



^ o o o o oo 
P oo o o o oo 



ro o 00 
o o o o 



ro 

o 







This figure "figl-4.png" is available in "png" format from: 


http://arXiv.org/ps/astro-ph/9502035vl 



COBE-DMR-NORMALISATION FOR COSMOLOGICAL CONSTANT DOMINATED 

COLD DARK MATTER MODELS 


Radoslaw Stompor^’^, Krzysztof M. Gorski^’^, and Anthony J. Banday'^ 


ABSTRACT 

The two-year CORE-DMR 53 and 90 GHz sky maps are used to determine the nor¬ 
malization of inflationary, cosmological constant dominated, cold dark matter models. 
Although the CORE-DMR data prefer a nearly critical value for the total density, low- 
density models with Oq ^ 0.15 can not be rejected at a confidence level 95%. We 
compare theoretical predictions for various large scale statistics, computed for a range of 
values of Oq and h = 0.5, or 0.8, to various measures of structure in the universe. 
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1. Introduction 


Cosmological constant dominated, spatially flat, cold dark matter models (CDM-A) 
can be considered a natural extension of the standard cold dark matter (CDM) scenario. 
By adjusting the combination of a non-zero cosmological constant and the total matter den¬ 
sity parameter, Oq (which remains poorly constrained by observations), we can maintain 
the spatial flatness of the universe required by the inflationary paradigm, whilst striving 
to improve the agreement of the theoretical model predictions to the observational data 
(Peebles 1984; Efstathiou, Maddox, & Sutherland 1990; Carrol, Press, & Turner 1992; 
Kofman, Gnedin, & Bahcall 1993). 

In this Letter^ we normalise the CDM-A models to the two year COBE-DM.Ii 
anisotropy measurements (Bennett et ah 1994), using the linear angular power spectrum 
estimation technique of Gorski (1994)f. Subsequently, we discuss the predictions of these 
appropriately normalised models for various large-scale structure measures. 

2. CMB Anisotropy Normalisation Procedure 

The data selection and power spectrum inference technique were implemented exactly 
as described in Gorski et ah (1994), and Gorski, Stompor, & Banday (1995). 

2.1 Theoretical Spectra of Anisotropy 

We consider spatially flat models with non-zero cosmological constant (Oq + A = 1) 
and primordial perturbations as prescribed (to within an arbitrary amplitude) by the 
inflationary scenario. Hence, random-phase, Gaussian, scalar curvature perturbations are 
assumed (no gravity waves), corresponding to adiabatic density perturbations with the 
Harrison-Zel’dovich power spectrum, P{k) oc k. The cosmological models considered are 
specified as follows: 1) The value of the Hubble constant. Ho = 100km s“^ Mpc“^, is taken 
to be given by either h = 0.5 or 0.8. 2) In agreement with Big Bang nucleosynthesis (BBN) 
constraints — = 0.013 (Reeves 1994) — we consider a baryon density O 5 = 0.05 for 

/i=0.5, or O 5 = 0.02 for h = 0.8. For comparison we also consider models with O 5 = 0.03 
for both values of h. 3) The total density, Oq, is varied between O 5 and unity and the 
value of A adjusted to retain a flat global geometry. 

The GMB anisotropy multipole coefficients and the matter perturbation transfer func¬ 
tions were computed using numerical Boltzman code as described in Stompor (1994) by 

f See also Bunn & Sugiyama (1994). 


2 



integrating the evolution equations until a redshift z = 0. Thus all effects caused by the 
cosmological constant were taken into account exactly. 

Figure 1 shows several of the CDM-A radiation power spectra normalised to the 
two year COBE-DM.Ii anisotropy. At low-t' ( 10), the radiation power spectra are de¬ 

termined by both the usual Sachs & Wolfe (1967, SW) and the cosmological constant 
induced, integrated SW (ISW) effects (Kofman & Starobinskii 1985, Gorski, Silk, & Vit¬ 
torio 1992), and depend very weakly on the Hubble constant. Over the range IQ ^ ^ ^ 30, 
an additional ISW effect, due to the non-negligible radiation energy density contribution 
to the total density after decoupling, is a significant factor in determining the shape of 
the anisotropy spectrum and introduces a noticeable dependence on h. The strength of 
both ISW effects falls with increasing Oq f- Consequently the power spectra flatten as Oq 
increases, and for models with near critical matter density the spectra monotonically grow 
with 7 . 

We express the power spectrum amplitude in terms of Qrms-PS ^ the value of the 
exactly computed quadrupole, 02 , multiplied by (S/dvr)^^^, which depends non-trivially 
on cosmological parameters. This is a generalisation of the normalisation introduced in 
Smoot et al. (1992) for pure power law model spectra. For Oq ~ 0.5, the ISW contributions 
boost the quadrupole over its pure SW value. For Oq ~ 0.5 — models with a negligible 
cosmological constant induced ISW effect — the resulting exact quadrupole is lower than 
its pure SW counterpart (see Table 1). The steepest spectra arise in models near critical 
matter density (Oq — 0 .8). 

2.3 Results of Qrms-PS Fitting 

Likelihood contours in the Qrms-PS — Oq plane (under the assumption of a uniform 
prior) are shown in Figs. 2a, b. The preferred value of the matter density is clearly close 
to critical, but even low-density models with Oq ~ 0.1 he well within the 95% confidence 
limits and appear acceptable. 

A typical likelihood fit for fixed Oq yields a ~ 13cr determination of Qrms-Ps{^ 
20 — 26/iK). Systematic effects, as discussed in Gorski, Stompor, & Banday (1995) and 

f However, in critical density models, only the cosmological constant induced ISW dis¬ 
appears. Furthermore, the ISW effect due to the radiation density contribution though 
usually small, is not negligible, and results in the steepening of the low-7 end of the power 
spectrum (Bond 1993; Bunn, Scott, & White 1994; Gorski, Stompor, & Banday 1995). 
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associated with 1) different coordinate frames (galactic/ecliptic) used for the (7(9i?-E-DMR 
data representation and 2) inclusion/exclusion of the quadrupole, are observed to shift the 
value of the best-fit normalisation by ~ 0.8 /xK and 0.4 /xK, respectively. Nevertheless, 
they do not influence any conclusions regarding the acceptability of the Iow-Oq models, 
and the lower limit of Oq ^ 0.15 is the most restrictive one. (Interestingly, this limit does 
not depend on h though the (7(9i?£l-DMR normalised Qrms-PS does.) 

A convenient summary of the inferred normalisation encompassing statistical er¬ 
rors, systematic shifts, and cosmological parameter dependences is then, Qrms-PS — 
{[20 -f 9.1/i'^ ® exp(— OQ/i'^'^/0.057)] x (1 ± 0.077) ± 0.4 ±0.2 ± 0.25}/iK. Exact values for 
a representative set of CDM-A models can be found in Table 1, which can be easily inter¬ 
polated to other models with high accuracy. However, our fit is good to within 5%. 

Note that, unlike in power-law or critical density models (Gorski et al. 1994; Gorski 
Stompor, & Banday 1995), there is no pivot point at a particular multipole order com¬ 
mon to all of the models considered. Nevertheless, for those models with ~ 0.05, a 

good approximate pivot can be found at / ~ 8 (Fig. 1), yielding the convenient model- 
independent normalisation og — 9.5 ± 0.73/iK. 

Since the dependence of the fitted Qrms-PS amplitude on model parameters is rather 
weak, and because the SW and ISW effects scale with the total density, there is a strong 
resulting variation with Oq of the amplitudes of the matter power spectra, as shown in 
Fig. 4. 

3. Measures of Large-Scale Structure 

Having determined the COBE-DM.Ii normalisations for our grid of GDM-A models, 
we now discuss the values of several large-scale structure statistics computed from the 
matter perturbation spectra. Hereafter, we apply the usual definitions: {'^hR)mass — 
(l/27r2)/“ dkPW]^HikR)P{k), and JsihR) = {R^/27r^)j\^ dkPWTH{kR) P{k), where 
P{k) is the matter power spectrum and Wth{x) = 3ji(x)/a;. Representative values of 
large scale structure measures derived from the models are included in Table 1. 

3.1 Mass Fluctuations: erg and Jg 

The high values of {(T%)mass > 1 predicted from the GORF-DMR normalised Oq = 
1, GDM models (Bunn, Scott, & White 1994, Gorski, Stompor, & Banday 1995) are 
difficult to reconcile with observations of the fluctuations in the number density of galaxies, 
{o's)gai — 0.8 — 1 (e.g. Davis & Peebles 1983; Fisher et al. 1993; Feldman, Kaiser & Peacock 
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1994; and Baugh & Efstathiou, 1993). At face value, this would require that the biasing 
parameter, b = {(T%)gaiI{(7%)mass ^ be somewhat smaller than unity, in direct disagreement 
with the estimates inferred from galaxy cluster catalogs (Henry & Arnaud 1991; White, 
Efstahiou, & Erenk 1993) and from the galaxy pair-wise velocities (Davis eA al. 1985, Cen 
& Ostriker 1992). A non-vanishing cosmological constant provides a means to circumvent 
these problems. The predicted (o's)mass as a function of is shown in Eig. 3a. A low 

value of Oo suppresses peculiar motions so high bias is not required to fit the small-scale 
galaxy pair-wise velocities. Thus unbiased models with (o's)mass — {c’'8)gal are likely to 
be of interest. However, physical motivations for requiring a biasing parameter somewhat 
larger than unity were discussed early on (e.g. Silk 1977), and are supported by recent 
hydrodynamical simulations (Katz, Hernquist & Weinberg 1992, Cen & Ostriker 1992) 
suggesting a value 6 ~ 1.5. In what follows, we conservatively discuss all values of 6 > 1. 
Models with ~ 0.12 then seem to be particularly interesting. 

Comparison of the Js-integral, computed on the larger scale of R = 20/i“^Mpc (and, 
therefore, more adequately modelled than (o's)mass by our linear approach), to the ob¬ 
servational value from Davis & Peebles (1983), J3(20/i“^Mpc) ~ 700 h~^ Mpc^ (probably 
accurate to within ~ 30%), also seems to favour the Iow-Oq CDM-A -models (Eig. 3b) in 
agreement with above constraint. 

3.2 Large-Scale Elows 

The rms amplitudes of large-scale bulk flows (Eig. 3c, d) are modelled as cr^(i7) = 
(27r^)“^i7g/(Oo) J dkP(k)W^p(k^ R)^ where Wbf = W^jj(kR)exp( — k‘^r‘^) is a model 
independent, POTENT-like window function as shown in Eig. 4 (r^ = 12/i“^Mpc). The 
factor /(Oo) ~ ® -|- (1 -|- Oo/2)A/70 as computed by Lahav et al. (1991) approximately 

describes the suppression of the velocity held in Iow-Oq CDM-A models). In Eig. 3d, the 
POTENT (Dekel 1994) and Lauer & Postman (1994) data are also shown. It is apparent 
that all models with % 0.1 can account for the POTENT velocities. However, if the 
Lauer & Postman data are confirmed even the critical density models will be in serious 
trouble. 

f In our calculations we smooth the velocity held obtained from the numerical integra¬ 
tion of the propagation equations, and these are the values given in Table 1. The Lahav 
et al. approximation is used for illustrative purposes, although it works well for Oq % 0 .1. 
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3.3 Galaxy and Matter Distribution Power Spectra 

Comparison of the theoretical matter power spectra with those inferred from the 
empirical data should, in principle, be the most complete and informative method possible. 
However, such a procedure is hampered by a number of factors, including the linear nature 
of our theoretical computations and the lack of agreement between different observational 
estimates of the galaxy power spectrum (which may be explicable in terms of intrinsic 
variations in the survey samples and/or by differences in the estimation technique utilised). 
In spite of such difficulties, we still consider such a comparison, shown in Fig. 4, to be 
informative in a qualitative way. 

On small scales {k > 0.1/i“^Mpc), the estimated power spectra from Baugh & Ef- 
stathiou (1994), Peacock & Dodds (1994) and the IRAS 1.2 Jy redshift survey are in 
reasonable agreement for those models with ~ 0.075. This limit is somewhat more 

restrictive than that obtained previously, and, combined with the COBE-DM.Ii plus PO¬ 
TENT inferred lower limits on Oq, leaves little room for a successful model. However, the 
results of da Costa et al. (1994) show enhanced power on scales k ki 0.1/i“^Mpc“^, con¬ 
trary to the expectation that such redshift space derived quantities underestimate the real 
space power spectra (Gramann et al. 1993). The observed discrepancies, particularly with 
respect to the similarly redshift-space-derived IRAS values, can not be solely ascribed to 
redshift space induced distortions. The relaxation of the above constraint would therefore 
seem to be in order until the observational situation is clarified. Conservatively, we retain 
the previously established relations. 

However, it then appears that the preferred models can not simultaneously fit the 
large-scale behaviour indicated by the APM data. This is more properly described by an 
unbiased CDM h = 0.5 model, which, unfortunately, fails noticeably on smaller scales. One 
should note, however, that on these larger scales there is again pronounced disagreement 
between different data sets (with the SSRS2-|-CfA2 data of da Costa et al. , 1994, clearly 
overshooting all other estimates by a factor ~ four, which cannot be completely explained 
by the effect of projecting the power spectrum to redshift spacef). It undoubtedly remains 

f This would require b ^ 20q'^®/ 5 < 1. Eor values of b in the range ~ 1 — 2, the linear 
correction appropriate for scales k ^ 0.03/iMpc“^ (Gramann, Cen & Bahcall 1993) is only 
10%-30% for Oo = 0.1 and 20%-50% for Oq = 0.4 so the pronounced power enhancement 
survives. 
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too early for any definitive conclusions to be reached. If the SSRS2+CfA2 data are con¬ 
firmed, the power enhancement visible on scales k ~ 0.01/iMpc“^ may favour the Iow-Oq 
CDM-A models. Improved, very deep galaxy catalogs may provide us with a definitive test 
of CDM-A models. 


4. Discussion 

In this Letter^ we reconsider the CDM-A models normalised to the two year COBE- 
DMR anisotropy data. We show that the Iow-Oq CDM-A (with Oq ^ 0.15) models can not 
be rejected solely on the basis of the COBE-DM.Ii data, although the most likely Oq value is 
strikingly close to unity (see also Bunn & Sugiyama 1994). Conversely, observations of the 
matter distribution do require low values of the total density resulting in the constraint 
^ 0.12. On larger scales, these models predict bulk flows in reasonable agreement 
with POTENT for ^ 0.1, but dramatically smaller than determined by the Lauer & 
Postman analysis for any (< 1) value of Oq- Therefore, if the latter observation is confirmed 
all CDM-A models (together with critical density CDM or MDM models) will be found 
wanting. The predicted mass distribution on scales of a few hundred Mpc overshoots 
most of the experimental estimates (including those derived from the IRAS 1.2Jy & APM 
surveys), but is in rough agreement with that recently obtained by da Costa et al. (1994). 

At present, although the position of the CDM-A models is more comfortable than that 
of critical density models, it is not free from potentially fatal flaws. Unfortunately, CMB 
anisotropies on a one degree scale do not offer any serious prospects for distinguishing be¬ 
tween cosmological constant models and other viable scenarios (Bond et al. 1994). Hence, 
definite conclusions will have to await more reliable observational data, particularly large- 
scale, deep galaxy surveys capable of unraveling the shape of the galaxy power spectrum 
on scales down to A; ~ 0.01 hM.pc~^ (ie. comoving length-scales up to 600 h ^Mpc). 

We acknowledge the efforts of those contributing to the COBE-DM.Ii. This work was 
supported in part by the Office of Space Sciences of NASA Headquarters. R. Stompor 
acknowledges LASP/COBE support during his stay at Goddard when this work was con¬ 
ducted. We are grateful to L. da Costa, K. Eisher and M. Vogeley for providing their 
power spectra, and A. Dekel for providing the bulk flow data. 
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FIGURE CAPTIONS 


Fig. 1.— CDM-A CMB anisotropy power spectra normalised to GORU-DMR. The thin 
lines correspond to the h = 0.5, O 5 = 0.05 models, thick lines to models with h = 0.8, 
O 5 = 0.02. The curves are drawn for four values of total density (in descending order at 
I = 2): Oo = 0.1, 0.2, 0.3,1. The approximate pivot point (for ~ 0.05 models) is also 
shown. 

Fig. 2.— a) Contour plots of the likelihood density in the Qrms-PS — plane for the 
h = 0.8 CDM-A models determined from analysis of the two year COBE-DM.Ii data 
in ecliptic coordinates. 68%, 95%, and 99.8% integrated likelihood levels, derived under 
the assumption of a uniform {Qrms-PS ^ Oo)-prior, and the Oq =const-ridge lines are 
shown. Thick lines correspond to the results of the analysis with the quadrupole included, 
thin lines excluding the quadrupole. Diamonds show the loci of the likelihood maxima: 
filled symbols show the situation when the quadrupole is included, open symbols when 
it is excluded. The filled and open circles show the likelihood maxima loci derived from 
analysis of the galactic data to illustrate the small shift in the inferred model amplitude. 
b) As in a for h = 0.5. c) Marginal likelihood densities for Oq, derived by integration of 
the two dimensional density distributions (from panels a and h) over Qrms-PS ^ are shown 
by thin lines. The functions corresponding to the analysis with the quadrupole included 
fall more steeply with decreasing Oq than when the quadrupole is excluded. The lines 
for the h = 0.8 model extend down to Oq = 0.02. Results derived with both the ecliptic 
and galactic frame data are shown — the ecliptic frame results fall below the galactic 
ones at small Oq- The heavy lines show the cumulative likelihood functions derived from 
the above likelihood densities. The quadrupole included case renders more steeply rising 
curves with h = 0.5, 0.8 models, and galactic/ecliptic results all merge together. The less 
steeply rising curves, split at very small values of Oq, correspond to the no quadrupole 
case. The 5% cumulative likelihood level is shown, as is the overall most restrictive, h- 
and coordinate frame-independent, Oq % 0.15 limit on the CDM-A models implied by the 
two year COBE-DMIi data alone. 

Fig. 3.— a) ((Ts)mass values predicted from the CORF-DMR normalised, CDM-A models, 
plotted as a function of . The h = 0.5, O 5 = 0.05 models are shown by solid lines 
terminating at = 0.25, while the h = 0.8, O 5 = 0.02 lines extend to the edge 
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of the plot. For each model, the upper of the pair of lines corresponds to the ecliptic 
frame/no quadrupole analysis, whilst the lower represents the galactic frame/quadrupole 
included case (the remaining two combinations would he in between). Dotted lines show 
the O 5 = 0.03 results averaged over the four possible combinations, b) The values for the 
J 3 integral over the density perturbation correlation function within 20 h~^M.pc (units 
[/i“^Mpc]^). Same coding as a. c) rms amplitudes of the bulk flow ([kms“^]) within 
spheres of 40, 60, and 100 h~^M.pc plotted in descending order, respectively, as a function 
of Qoh. Heavy lines — h = 0.8, Oq = 0.02, thin lines — h = 0.5, Oq = 0.05. d) rms 
amplitudes of the large scale flows ([kms“^]) as a function of the top-hat sphere radius 
encompassing the sample (an additional = 1200 km s“^ smoothing is included). Heavy 
lines — h = 0.8, O 5 = 0.02, and Oq = 1, 0.4, 0.3, 0.2, 0.1, 0.05, and 0.02 in descending 
order. Thin lines — h = 0.5, O 5 = 0.05, and Oq = 1, 0.4, 0.3, 0.2, 0.1, and 0.05 in 
descending order. The POTENT data (courtesy A. Dekel) are shown by circles, and the 
square denotes the Lauer & Postman datum. 

Fig. 4.— COHE-DMR normalised inhomogeneity power spectra. The h = 0.5, O 5 = 0.05, 
and h = 0.8, O 5 = 0.02 models are shown in thin and heavy lines, respectively. Dotted 
lines correspond to Oq = 1, CDM models. Solid lines, in descending order, show the 
CDM-A power spectra for Oq = 0.1, 0.2, 0.3 — heavy type, and Oq = 0.1, 0.2, 0.3, and 0.4 
— light type. Overplotted are several observational estimates for the galaxy distribution 
power spectrum: squares - Baugh & Efstathiou (1993), filled circles - Peacock & Dodds 
(1994), pentagons - Fisher et ah (1993), triangles - da Costa et ah (1994) (redshift 
space). Miscellaneous spectral windows required for the computation of the statistics 
considered in this paper are shown (as annotated) as the lowermost curves. The quadrupole 
window depends (due to the SW and ISW effects) on model parameters, Q'^ms-PS ~ 
(27r^)“^ j dk P{k)WQ{k). For the pure Sachs-Wolfe effect, WQ{k) = j|(A:i? 2 ^ 55 )/ 4 A:^, 

where Rlss is the radius of the last scattering surface (at a redshift zj^ss — HOO). The 
windows for Oq =0.1, 0.3, and 1 are shown from left to right (note the trend in the window 
maxima and changes in the shape of the functions). 
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